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ABSTRACT 


A mulitary aircraft deployment problem from the United States Transportation 
Command is modeled as a generalized transportation problem with side constraints and 
solved using a general purpose linear programming package. The model involves the 
assignment of military units and material to aircraft and the assignment of aircraft to 
missions in order to appraise the utility and to determine the assets required for prelim- 
Inary military operation plans. A transformation of this model which aggregates vari- 
ables relating to the early or late delivery of requirements is also described. A specialized 
algorithm which separates an instance of the model into subgroups of independent time 
windows, finding the globally optimal solution by solving independent subproblems 1s 
also explored. Lastly, an integer rounding model is described which converts continuous 
solutions to integer in order to facilitate implementation of the former models with an 
existing post-solution processor. 

Excellent quality solutions are provided for problems involving nine routes, 80 
movement requirements distributed across two cargo classes involving 200.000 short 
tons of freight. and 250 aircraft using four different aircraft tvpes for each of I2 time 
periods. The problem. which has the potential of having over 10,000 variables. ts re- 
duced significantly using variable reduction and the aggregation transformation. The 
reduced problem requires approximately 1,000 variables and 300 constraints and sol- 
utions are obtainable in under 14 minutes using the General Algebraic Modeling Svstem 


on an 80286-based personal computer. 
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I. INTRODUCTION 


In the modern world of rapid mobility, the difference between a successful and an 
unsuccessful deployment of miltary units and equipment will depend on the efficient 
allocation of transportation assets to the requirements which need to be moved within 
prescribed time intervals. These transportation assets consist primarily of mulitary 
transport aircraft, ocean-going vessels, and rail and truck vehicles. The movement re- 
quirements may consist of any material item or person requiring transportation from one 
location to another such as an infantry division, petroleum products, heavy equipment, 
artillery, subsistence items, and many others. Small differences in the timing of delivery 
or the effective use of assets could deterniune the difference between an action being won 
omiost |Ref. I}. 

This thesis considers the question of whether a potential military operation can be 
Initiated and sustained from the perspective of air transportation. It 1s valuable to know 
What assets will be used, where thev will be needed and whether there are bottlenecks in 
the transportation process which significantly affect the operation. It 1s not necessary 
to develop detailed schedules and load plans for aircraft and movement requirements for 
this mitial planning: creating such schedules and plans 1s bevond the scope of this work. 
However, the initial planning process for the allocation of transportation assets will need 
to be sensitive to the timing of delivery of these movement requirements and to the ef- 
fective use of the transportation assets. The problem of determining the appropriate 
deplovment of aircraft given a set of requirements Which must be moved within pre- 
scribed time intervals and with specific delivery umes which are the most opportune for 
the objectives of the operation will henceforth be called the deployment problem. Past 
model development has shown that the complenities involved with combining different 
modes of transportation in a single model require solving an extremely large problem. 
Consequently. a modular approach which separates the air, land. and sea allocation 
problems is usually emploved [Ref. 2]. 

The United States Transportation Command (USTRANSCOM), is the Department 
of Defense’s agency which has been most concerned with the deplovment problem. This 
organization was created to consolidate the global air, land, and sea transportation as- 
sets and requirements which are used by the military to meet national security objectives. 


While tasked with coordinating the development of deployment plans tmvolving 


transportation assets, one of USTRANSCOM(/s principal endeavors 1s to provide 
automated data processing svstems to streamline the planning process. The Joint 
Deplovment Agency (JDA), the predecessor organization of ULSTRANSCOM, 
recognized the potential of automated data processing to assist with mobilization plan- 
ning. The Joint Deployment System (JDS), a comprehensive computer planning system 
initiated under the auspices of the JDA, was designed to plan, execute, and monitor the 
deployment of units during both peacetime and times of crisis. Manv reseach initiatives 
directly related to the deployment problem have been undertaken to support this system. 
[Rete 

The first linear programming based effort at modeling and solving the deployment 
problem was the System for Closure Optimization Planning and Evaluation (SCOPE), 
a submodule of the JDS [Ref. 4]. This system used a Benders decomposition method to 
solve a continuous Variable model consisting of both sea and air assets over discrete time 
periods. The formulation separated the model into a master problem involving the as- 
signment of assets to delivery missions and a subproblem involving the assignment of 
movement requirements to transportation assets. Because the solution method con- 
verged slowly or not at all. and because it did not model sea assets integrally, the Benders 
decomposition approach was considered to be unsuitable for the JDS svstem [Ref. 5: 
De 20) 

The deployment problem for the allocation of sea assets has been studied by manv 
researchers since the original SCOPE model. A continuous variable deplovment model 
for both sea and air assets similar to that of the SCOPE model was undertaken by 
Colher [Ref. 4] using hnear programming and a variable reduction heuristic. This model 
was developed to assist planners with designing feasible deployment plans and selecting 
appropriate modes of transportation for movement requirements. An alternative model 
for the deplovment problem involving only sea transportation assets was developed by 
Lally [Ref. 6] using integer programming. This model avoided the integrality problems 
Which occurred for sea assets in the SCOPE model. While both of these studies dem- 
onstrated the ability to use optimization techniques for moderately sized instances of 
problems, neither of these techmques have been employed in the JDS system. In its 
current form. the JDS uses simulation-based heuristics which are capable of providing 
solutions to highly detailed and very large resource allocation models. However, the 
JDS still has difficulty in providing solutions to large instances of problems in a timelv 
manner [Ref. 7], and the solutions obtained cannot be considered to be optimal with 


respect to any measure of effectiveness. More promising developments involving opti- 


to 


mization techniques for solving large-scale deployment problems involving only sea 
transportation assets have been made by Buvik [Ref. 8] and Lima [Ref. 9] using set 
covering models with dynamic column generation. 

Another initiative Which is being pursued by USTRANSCOM to solve the deploy- 
ment problem is the development of the Deployment Analysis Prototype (DAP). The 
DAP is a deployment modeling system developed principally as a technology demon- 
stration implementable entirely on a $80286-based personal computer using manv of the 
same inputs as the JDS. This svstem is used as a base for the development of modeling 
techniques which could be used as stand-alone tools for deployment planning as well as 
a testbed for the development of innovative planning techniques for the JDS. Of specific 
interest for this thesis will be the continuous variable model designed to solve the de- 
ployment problem for air transportation assets which was developed by Rosenthal 
[Ref. 10] to succeed the original air transportation allocation model in the DAP. 

The nature of deplovment modeling on a large scale poses many challenges when 
optinuzation methodology is used. This thesis develops prototypic planning tools to aid 
USTRANSCOM with the allocation of air transportation assets in an effective and 
timely manner. Formulations are developed using the same problem format. but with 
a different modeling approach, of those previously laid out by Rosenthal for the allo- 
cation of air assets. The purpose for developing these new models 1s to enhance system 
solution-time performance, model resolution, flexibihty, and compatibility with other 
modeling programs. Three optimization models and an iterative solution algorithm de- 
signed to improve the DAP are provided. The models and algorithms form a collection 
of ideas related to the deployment problem as it relates to air transportation. 

Chapter 2 of this thesis provides a general discussion of the aircraft deployment 
problem as Well as the derivation and description of the optimization models and algo- 
rithm. The first model is based on the weighted distribution problem, more commonly 
known as the generalized transportation problem [Ref. 11]. As in earlier models, this 
model uses continuous variables to represent the assignment of movement requirements 
to aircraft and the assignment of aircraft to missions. The convention of using contin- 
uous variables is justified by the continuous nature of air operations, the aggregate na- 
ture of available data and the level of detail required for the preliminary screening of 
potential operation plans. Unlike the earlier models, this model provides information 
on the movement requirements which are not deliverable at the opportune time for the 
operation plan being considered. This information is particularly important for deciding 


Whether sufficient deliveries of particular movement requirements will occur and thus 


allow the military operation to proceed. Additionally, this model lays the foundation 
for developing the other models. The second model, a transformation of the generalized 
transportation formulation, significantly reduces the number of variables when move- 
ment requirements exist across most route, opportune delivery period, and cargo type 
combinations. This model exhibits a more stable model generation and solution struc- 
ture than the generalized transportation formulation when solved as a linear program. 

A solution algorithm which searches for and then solves subproblems over inde- 
pendent time windows of the original time horizon is described. In many instances, this 
technique may offer the ability to model extended time horizons and solve larger in- 
stances of problems. 

The DAP, in its current form, uses a post-solution processor which requires integral 
aircraft solutions. In order to reform the linear programming solution to integers as 
required by this post-solution processor, the resulting nonintegral solutions of the con- 
tinuous variable models are converted to integer solutions. This third model uses a 
greedy rounding model to allocate fractional values making use of any remaining aircraft 
supplies. 

Test results of typical problems are provided in Chapter 3. A typical problem size 
involves nine routes. 80 movement requirements distributed across two cargo classes 
involving 200,000 short tons of freight, 250 aircraft per time period using four different 
aircraft types, with 12 time periods in the time horizon. A summary of the thesis 1s 
provided in Chapter 4. All of the models and techniques described in this thesis are 
programmed in the General Algebraic Modeling Svstem (GAMS) [Ref. 12] using a 80286 
based personal computer. GAMS was used because it facilitates rapid changes in the 
model constraints to reflect unrque situations for a particular nulitarv operation. 


GAMS programs for the models described in Chapter 2 are contained in the appendices. 


II. PROBLEM FORMULATION 


A. GENERAL PROBLEM DESCRIPTION 

The military aircraft deplovment problem involves the shipment of military units 
and equipment between ports of embarkation and debarkation within prescribed time 
intervals using an available pool of aircraft. These shipments are assumed to be 
unidirectional, going from one geographic area to another. An example of this would 
be the movement of a rapid deployment force from the east coast of the contimental 
United States to Europe. This chapter presents the derivation of a collection of models 
and an algorithm which can be used to evaluate the general suitability of preliminary 
military operation plans which depend on the solution of this problem. The remainder 
of this sectron provides a general description of the terminology and elements of the 
models. 

For the models used in this study, the planning horizon 1s a collection of discrete 
non-overlapping time intervals called periods. The mnterval of a period depends on the 
distances to be travelled and the speed of the aircraft. Any cargo shipped in a period 
will arrive at its destmation within that same period. Additionally, the mtervals for the 
different periods may be of equal or varied duration. One reason for allowing the du- 
ration of time periods to vary is because of the increasmg uncertainty of events and re- 
quirements the farther the trme horizon 1s extended into the future. Consequently. the 
amount of detail required in the mode! dinunishes with time which allows increasing the 
trme interval as periods progress. 

The deployment problem involves the transfer of military units and equipment. Any 
collective group of materials or tndividuals with a common port of embarkation and 
debarkation which is requested to be moved within a specific time period is called a 
movement requirement. The period of requested delivery is called an opportune delivery 
period since it reflects the period of delivery which is deemed to be most beneficial for 
the objectives of the operation plan being considered. In addition to the opportune de- 
livery period, an interval of permissable periods is defined which contains those periods 
for which movement requirements are allowed to be delivered early or late from the 
Opportune delrverv period. An available-to-load period is also used which indicates the 
earliest possible period for which a particular movement requirement can be loaded on 


ai aircraft. This period represents the availability of cargo, while the earlv delivery 


period mentioned above represents the earliest delivery period allowable which still 
meets the objectives of the operation. 

All items within a particular movement requirement will have a single homogeneous 
classification with a standard physical unit. For clarity, the definition of the term phys- 
ical units is noted as being quite distinct from the term military units. While the term 
military units refers to an organized bodv of troops, physical units refers to the fixed 
quantity of measure for items being transferred. The classifications of movement re- 
quirements are typicallv passengers, bulk cargo, oversized cargo, or any other appropri- 
ate grouping, while the physical units of measure are typically tons for equipment and 
passengers for mulitary units. These movement requirement classes will be called cargo 
for simplicity. 

The aircraft which can be deployed during a particular period are not treated indi- 
vidually, but rather as a pool of aircraft available to perform missions. The assumption 
that deplovable assets are pooled without regard to the location of aircraft at the start 
of a period is reasonable because of the continuous nature of aircraft movements and 
the assumption that all embarkation ports are in the same geographic region. This 
simphfication significantly reduces the number of variables required in formulating the 
model since the location of each aircraft does not have to be tracked. Anv of the aircraft 
assigned to the pool mav be used on successive tours between a single port of 
embarkation and debarkation during that same period. These tours are designated by 
the forward routes which connect the port of embarkation with the port of debarkation. 
Aircraft missions are defined to be a route, delivery period, and aircraft tvpe combination 
which are used to transport various movement requirements. For this study, the aircraft 
deliverv periods are called mission periods. Since the deplovment problem assumes a 
unidirectional flow of military units and equipment. aircraft are considered to return 
empty from ports of debarkation. Aircraft are classified bv aircraft-twpe and have asso- 
ciated movement requirements which are compatible with the aircraft design. The air- 
craft are assigned from the pool based on the relative importance of the aircraft. One 
reason for assigning aircraft usage priorities is that certain aircraft are better suited to 
particular military operations. 

The aircraft deplovment models of this study use formulations with a channel flow 
concept similar to that originally proposed in the SCOPE model [Ref. 3]. Channels or 
routes are defined only on possible pairs of embarkation and debarkation ports where 
movement requirements exists. The capacity to transport quantities of a requirement 


over a route using a specific aircraft tvpe is calculable. This route capacity is a function 


of manv factors such as the expected aircraft speed, the length of the interval of time in 
the period, the distance between ports of embarkation and debarkation, the loading ca- 
pacity of the aircraft. the expected lavover time at ports, and the compatibility of aircraft 
tvpes with particular ports. The route capacity is an input parameter available to the 
model. 

The model formulations employ continuous variables to represent shipments of 
movement requirements across routes. A shipment assigned to a mission is defined to 
be delivered entirely within the period of the mission. The quantity of cargo shipped per 
period on a route using a particular aircraft is determined by its route capacity. As such, 
shipments are based on the rate of transfer specified by the route capacity during a pe- 
riod and not discrete plane-loads of cargo. The amount of cargo transferred is limited 
by the total supply of aircraft available. The use of continuous variables is reasonable 
when nontrivial numbers of aircraft are required on any route during a military opera- 
tion. For such cases, the continuous solutions should be sufficiently accurate to provide 
a good preliminary assessment of the logistic supportability of an operation plan. The 
models are required to be able to not only track the closure of all movement require- 
ments to their debarkation ports. but also those portions of movement requirements 
Which are delivered earlier or later than their opportune delivery period. This is espe- 
cially important for determining whether a marginal operation plan is supportable. 
Aircraft supply constraints are imposed on the model bv using total route capacity re- 
strictions determined by the total number of a particular aircraft tvpe available during 
a period. Various linear side constraints can be included to adapt the model to specific 
operational situations such as capacity limitations at ports or minimum shipment level 
restrictions for movement requirements. 

The objective function consists of costs which are combinations of penalties assigned 
to the different events associated with transporting movement requirements. The cost 
scheme for a model includes relative penalties for (a) the nondelivery of movement re- 
G@mmemients, (b) the early or late delivery of cargo from the opportune delivery period, 
and (c) the assignment of cargo and aircraft tvpes. The additive combination of the 
different classes of penalties creates a cost scheme which prioritizes all possible ship- 
ments which could occur using the various movement requirements, aircraft, and deliv- 
ery periods. The specific cost of a transfer may not be economically meaningful: 
however, the relative cost of the transfer in relation to the cost of other types of transfers 
is important. The dominant penalties typically are for delivery deviations from the 


opportune delivery period and nondelivery of movement requirements. The different 


tvpes of penalty classes are hierarchically scaled so the dominance of particular penalty 
classes 1s preserved when they are added to lower priority penalty classes. While the 
objective function minimizes the costs derived from these penalties, the general objective 
of the model ts to deliver all movement requirements mm a timely manner while efficiently 


using available aircraft. 


B. GENERALIZED TRANSPORTATION FORMULATION 

The original formulation by Rosenthal [Ref. 10] for the aircraft deployment problem 
used dynamic equations to model the shipment of movement requirements and early and 
late deliveries aggregated by mission delivery period in order to reduce the problem size. 
However, thts model was unable to provide information about the or1gination of move- 
ment requirements which were delivered early or late and did not consider the possibility 
of an infeasible problem. The models of this thesis take these aspects of the deployment 
problem mto consideration with emphasis berng place on the size and solvability of 
problems. This section presents a model which solves the aircraft deployment problem 
using a weighted distribution formulation [Ref. 11]. 

The werghted distribution formulation is a generalized transportation network which 
weights the variables in a problem with non-unit coefficients. For the aircraft deploy- 
ment problem. these weighting coefficients are multipliers that represent the route ca- 
pacities Which regulate the transfer rates of movement requirements. Side constraints 
may be added to the problem to represent specific operational situations. Many of these 
side constraints may be formulated so thev extend the network of this Generalized 


Transportation Formulation (GTF). The primary index sets for the model are: 


tel embarkation ports. 
jeJ debarkation ports. 
reRkR routes. 

aéea aircraft types. 
heH Careonty pes. 


i,7¢P ~ periods of time horizon: 


For notational purposes, ¢ and t are dummy indices representing any period in P. The 
location of ¢ and t in a data parameter or variable determines the context for the use of 
either of these dummy indices. The first index set involving P in a data parameter or 


variable always relates to a mussion period while the second index set involving P always 


relates to the opportune delivery period. Additionally, the set of routes R , represents 
the cartesian product of embarkation and debarkation ports. 

For simplicity, the data for the problem is expressed using a single measurement 
system in short-tons (stons), stons’aircraft, etc., which are used in a specific period. 
Measurement systems using volumetric or other physical units are possible. The data 


for the problem is: 


MIREQ_, movement requirement of cargo tvpe A in stons to be shipped over route r 
with opportune delivery period Tt . 


ALP_, available-to-load period for a movement requirement . 


meds, number of periods which a movement requirement may be delivered early 
or late from its opportune delivery period. 


ASUP,, number of aircraft of tvpe @ available in mission period 7. 


RCAP.,, cargo capacity tn stons, aircraft for aircraft type a@ over route r during 
mission period 7. 


PCAP,, — port throughput capacity for movement requirements involving cargo type h 
in stons for port / during mission period 7. 


ce. cost per ston for shipping cargo tvpe with opportune delivery period t 
over route r using aircraft tvpe a during mission period 7. 
Uh cost per ston for not delivering cargo within a movement requirement hav- 


Meacencoum PeyP Avi Opportiine deliver. perod < required Over route 7 . 


The satisfaction of movement requirements 1s accomplished using variables associ- 
ated with avatlable aircraft. The delivery of the movement requirements can be identi- 
fied as belonging to two different classes of assignment. The first class of assignment 
represents the direct satisfaction of a movement requirement with aircraft having a 
mission period the same as the opportune delivery period, henceforth to be called direct 
period deliveries. Ihe second class represents the satisfaction of movement requirements 
by aircraft with a different mission period than the opportune delivery period (vielding 
early or late deliveries of cargo) henceforth to be called cross-period deliveries. Both of 
these classes are represented in the variable X,,,,. If there are inadequate supplies of 
available aircraft to satisfy all movement requirements, variables associated with 
undelivered cargo are used. The variables for the undelivered movement requirements 
will be called infeasibility variables. These variables bear a sufficiently large cost so there 
use occurs only when there is no other means of delivery. The variables for the model 


ale: 


an stons of cargo type A with opportune delivery period t shipped over route r 
using aircraft tvpe @ during mission period 7. 


N 


rth 


undelivered stons of a movement requirement having cargo type 
A with opportune delivery period t required over route r. 


In addition to the primarv index sets, the following derived index sets are necessary. 
The set of periods P in the time horizon must contain the set of periods P,, for which 
aircraft are available to perform missions, that is P,< P. The set of periods P in the 
tume horizon must also contain the set of periods P, for the opportune delivery periods 
of the movenient requirenients, tiidtse ater 

The set Z of all combinations of routes, opportune delivery periods, and cargo types 


for movement requirements Is: 
LS eke ae a 


The set LCL of allowed combinations depends on the existence of a movement re- 
quirement involving a particular route, opportune delivery period, and cargo type com- 
bination. The set S of all combinations of aircraft tvpes available in a given mission 


period 1s: 
Sara x Pea 


The set Sc S of allowed combinations depends on whether there are any aircraft avail- 
able during a given mission period. 
The set Q of all possible combinations for transporting movement requirements on 


different aircraft types during mussion periods 1s: 
O = {(AxP, é€ S) ON Uo ete Dy, 


The set 0 c Os the allowed combinations of A.R, P, , and H such that nonzero 
route capacity exists and a compatible aircraft and cargo combination exist. Addi- 
tionally, the mussion period in P,, must fall within the interval containing the early devi- 


and late deviation d,,, from the opportune delivery period in P, and must also 


rth 


ation @., 
be later than the available to load period ALP,,, for a particular movement requirement. 
This set ensures the combination of primary indices occurs only when there is a move- 


ment requirement and the aircraft and route capacity exists to deliver the requirement 


10 


Within mussion periods which fall inside the interval of permissible periods. The set Q,, 


of all movement requirements using a particular aircraft type @ and mission period 1 is: 


ee 


O, = {Rx P,xH | Rx P,x Hx {(a)} € Q}. 


The set Q,,, of all aircraft tvpes available during the mission periods P,, with a particular 


route ry , opportune delivery period t , and cargo tvpe /iis: 
O ch = { A x Er | 7 Pn x CU aicenz)) € Q }. 


In addition to the conditions mentioned above, L, S , and Q may contain other implicit 
constraints which are specific to a given problem. The model is formulated as a hnear 


program: 


aa we? a : Noy 
Miuninuze ae Chat naa » Crinaan 


(Qy,f, -n)e O (r,7,h)e L 
Subject To 

De /RCAP,, = ASUP,, Ges dep 

(r, Te O,, 
ee Nee REO. Gaie L anee, 

(@,0)€ Om 
linear side constraints involving X,,,,-; (ies) 
O< Narich (1.4) 
rs (1.5) 


rth 


Equation (1.1) ensures that all movement requirements are satisfied bv actual slup- 
ments or infeasibilitv variables. Equation (1.2) ensures the number of aircraft deploved 
for missions does not exceed the available pool of aircraft. Equation (1.3) represents 
linear side constraints such as those associated with port throughput restrictions. An 


example of such constraints for ports of embarkation 1s: 


ve, < PCapaem (1.3’) 


Lartch 
(aca) 


(r=(Z/) | Ceti }xJ ) 


A similar set of constraints could be constructed for the debarkation ports. Equations 
(1.4) and (1.5) are standard nonnegativity constraints. 

The first two constraint sets, equations (1.1) and (1.2) form the generalized trans- 
portation network. It 1s noted that port capacity constraints (1.3’) could be included as 
generalized network constraints by expanding the transportation network to a 
transhipment network. This is not pursued here since not all types of side constraints 
can be handled 1n this manner. A simple three period depiction for the bipartite graph 
of the generalized transportation network 1s shown in Figure 1. Nodes are divided into 
those related to the supply of aircraft of a particular type available during a mission 
period (S) and those related to the movement requirements (L) . The arcs, represemteg 
by the variables, are directed from Sto L. Direct period delivery arcs are shown with 
solid lines while the cross period delivery arcs are shown with dashed lines. 

The cost scheme for the variables establishes a set of hierarchical penalties based 
On various aspects of the deployment problem. The primary tier of the cost scheme in- 
volves the timeliness of delivery for the three different classes of assisnment., “Vinememe 
no penalty for the direct delivery of a movement requirement. Cross period delimem 
penalties of PPEN, are assessed dependent upon whether a delivery is late or early and 
upon the amount of the deviation from the opportune deliverv period. The penalty 
structure for cross period variables is an increasing function of the earliness and lateness 
of movement requirement delivery. This function imposes a precedence pattern on the 
model which imphes a delivery of one day early 1s better than one dav late which 1s better 


than two davs early and so forth. The following relation illustrates this priciple: 


0 = PPEN., < PPEN 2 =< PPEN.,, . < PPEN.). 


This pattern progresses throughout the entire interval of permissible periods for delivery. 

Two other penalty structures subordinate to the first are imposed, which allow fur- 
ther refinements of the objectives for assigning aircraft to movement requirements. 
These two additional penalty structures form the second and third tiers of the penalty 
scheme and are based on the relative priority to be used for assigning cargo and aircraft 


tvpes to missions. The second tier of the penalty scheme assesses penalties of //PEN, 
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Figure J. Generalized Transportation Network 
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based on the assignment priorities for cargo classes involved with cross period deliveries. 
This penalty tier is used to prioritize cargo in the event of an early or late deliverv. A 
tertiarv tier of penalties APEN, , assesses penalties for using particular types of aircraft 
for any type of feasible deliverv. This penalty tier is used because certain types of air- 
craft may be considered to be more suitable to the military operation being considered, 
while other aircraft mav be more valuable as a reserve component. 

The hierarchical penalty structure gives a precedence relation for assessing the im- 
portance of variables when higher level penalties are equal. It 1s noted that each of the 
three penalty tiers are defined over index sets which are independent of the mndex sets 
of other penalty tiers. Therefore, the individual penalties within a penalty tier are spe- 
cifically related to its index set. Additionally, all subordinate penalties have a sufficiently 
low value so that they do not influence the penalty structure of higher tier penalties, that 
is, lower tier penalties are smaller than the smallest difference of penalty magnitude be- 
tween anv two penalties associated with a higher tier. This creates a cost scheme which 
allows the additive combination of priorities without entangling the various objectives 
represented by the precedence relations of the different tiers. The costs associated with 


particular variables is derived form combining penalties as follows: 


C.., = PPEN, + HPEN, + APEN,. 


ariczh 


Additionally, the cost for infeasible deliveries of Cx, 1s set larger than the maximum of 
the feasible variable delivery costs in order to ensure that infeasibilitv variables are used 


only if feasible delivery is impossible. 


C. MOVEMENT REQUIREMENT REALLOCATION MODEL 

The Generalized Transportation Formulation may have a large number of variables. 
A large problem for an 80286-based personal computer might involve twelve time peri- 
ods. four aircraft types. two cargo types. and nine routes. There 1s the possibility that 
over 10,000 variables could be generated which would most Irkely make the problem 
intractable to solve using a general linear programming package. For this reason, a 
transformation of the Generalized Transportation Formulation 1s proposed in order to 
reduce the number of variables. 

In the original formulation, variables are required to track each possible movement 
requirement that could be assigned to various missions. An alternative method for 
tracking assignments is to create separate variable types for the quantity of cargo 


amongst all opportune delivery pertods assigned to each mission and for the quantity 
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of movement requirements which are delivered in periods other than the opportune de- 
livery period regardless of aircraft tvpe. With this alternative method of variable gener- 
ation, the direct assignment of cross period deliveries of aircraft missions to movement 
requirements is avoided. This limits the size of the problem mentioned above to under 
3700 variables. This new model is equivalent to the Generalized Transportation For- 
mulation with following exception. In the Generalized Transportation Formulation, the 
tvpe of aircraft was determined for the cargo of movement requirements which was de- 
livered early and late as well as that which was delivered on time. In the new model, the 
cargo which 1s delivered early or late will be reallocated to cargo of some other period's 
movement requirement. Therefore, this late or early cargo 1s not specifically tracked as 
a separate delivery on a particular aircraft type, but rather as a part of a different 
movement requirement. The information lost 1s the aircraft tvpe which the early or late 
cargo arrived on. This information 1s usually irrelevant since it does not matter which 
@inevairt the late or early cargo arrived on, but rather which cargo armved early or late 
and how early or late it was. The transformation is unplemented using an aggregation 
of shipments across opportune delivery periods for movement requirements using a 


particular route and cargo type delivered during a specific mission period as follows: 


vs ae. Ne 
iy | af Fate? at ae 


f 


Additionally, aircraft tvpes are aggregated for each mission period involved with a 


movement requirement as follows: 


ee = Dy eh (J ig Oe 


Qa 


The transformation of X,,.., into these two new variable types is derived by separat- 
ing the possible variables into direct period deliveries and cross period deliveries. The 
direct period deliveries can further be separated into a positive term involving all possible 
cargo shipped on a particular mission and a negative term involving all early or late 
cargo shipped on the mission. Algebraically, the derivation of the transformation for a 
feed (7, t, Jt) 1S: 


» > Ager a » » er ae » y Spien 
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The final right hand side of this formula has three terms. The first term represents the 
aggregate movement requirements which are transported on a mission using aircraft tvpe 
ain period t . henceforth to be called the mission delivery variable. The second term 
represents the amount from each of the movement requirements not having an 
opportune deliverv period of t , Which are transported in mission period zt. This is a 
reallocation of movement requirements from periods other than t by regrouping them 
With a movement requirement which has opportune delivery period t. The third term 
represents the reallocation from the movement requirement with an opportune delivery 
period of t by regrouping them with movement requirements having opportune delivery 
periods other than t. The variables associated with these last two terms will be called 
the reallocation variables. 

When the infeasibility variable .V.., and the movement requirement ./REQ.., are 1n- 
cluded with this transformation, a restructured version of equation (1.1) of the Gener- 
alized Transportation Formulation is obtained. A conservation of flow relation exists 
for a specific movement requirement node involving the set of indices (r,t. A). Simi- 
larly, A”,,., 1s substituted into the remaining equations to give commensurate transformed 
versions of the original equations. The new formulation is called the Movement Re- 


quirement Reallocation Model (MRR) The primary index sets for the model are: 


ie] embarkation ports. 


jeJ debarkation ports. 
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reR routes. 
aeaA aircraft types. 
heH cargo types. 


fmer =§©6penods of time horizon. 


The data for this model uses the same physical units convention as the generalized 


transportation model. The given data for the problem 1s: 


A{REQ.,, movement requirement of cargo type A in stons to be shipped over route r 
With opportune delivery period t. 


ALP_, available-to-load period for a movement requirement . 


— 


d_.,a., number of periods which a movement requirement may be delivered early 
or late from its opportune delivery period. 


ASUP,, number of aircraft of tvpe a available in mission period ¢. 


RCAP,,, cargo capacity in stons’aircraft for aircraft type a over route r during 
mission period r. 


PCAP,, port throughput capacity for movement requirements involving cargo type /t 
in stons for port / during nussion periods. 


ee; cost per ston for shipping cargo type / with opportune delivery period 1 
over route r using aircraft tvpe @ during mission period r. 
Ge. cost per ston of reallocating a mission requirement of cargo type /t over 
foule we Trem a Opponiune deliver period of «to one of ¥ . 
Bs cost per ton for not delivering cargo within a movement requirement having 
cargo type A with opportune delivery period t required over router. 


The variables for the problem are: 


ACB stons of cargo type A shipped over route ry using aircraft type a during 
mussion period 7. 

ee, stons of a mission requirement with cargo type A over route r with 
opportune delivery period t which is reallocated to mission period 7. 

ia undelivered stons of a movement requirement having cargo type 


h with opportune delivery period t required over router. 


In addition to the primary index sets, the following derived index sets are necessary. 


The sets P_ , P,. and S are defined the same as in the General Transportation Formu- 


7, 


lation. The set L’ of the cartesian product of routes, periods, and cargo types for 


movement requirements 1s: 
C= 1 hee 


The set L’¢ L' contains the allowed combinations of route, period and cargo indices. 
At least one of two conditions must occur for an index combination to be a member of 
L’. Either a movement requirement involving a route, the opportune delivery period, 
and the cargo type exists, or there are aircraft types available to deliver movement re- 
quirements during the opportune delivery period. The last condition does not require a 
nonzero movement requirement to exist since the model allows the reallocation of com- 
patible movement requirements from other opportune delivery periods. 

The set }’ of all possible combinations of indices for transporting movement re- 


quirements on different aircraft tvpes during the possible mission periods 1s: 
=o BAe eee ae 


The set }'c V restricts this set by allowing only combinations a, r, re P, hk such that a 
nonzero route capacity exists and there are aircraft available which are compatible with 
the cargo tvpe. The set IF of all possible combinations of nussion periods with move- 


ment requirements 1s: 
Po a aX ee 


The set J’ c Ho of allowed combinations contains only indices which have nonzero route 
capacity for at least one of the available aircraft tvpes of a specified movement require- 
ment. The aircraft type must be compatible with the particular cargo type of #7. Ad- 


diuonally, the mission period in P,, must fall within the interval containing the early 


m 


deviation d., and late deviation 


srt); 


.., rom the opportune delivery period in P, and must 
be later than the available loading period ALP..,. This set ensures that the combination 
of primary indices occurs only if a movement requirement exists and there are compat- 
ible aircraft and route capacity available which could deliver it using mission periods 
Which fall within the interval of permissible periods. 

The set ae of all aircraft tvpes which can be used for a particular movement re- 


quirement 1s: 


ne ee! 


Vv, = {A | Ax f(r, 7h)} € V3. 
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The set V,, of all routes and cargo tvpes which can be used with a particular aircraft tvpe 
and mission period combination is: 


{Rx H| Rx Hx {(a)}e V}. 
The set H’_, of all mission periods for a particular movement requirement 1s: 
ae een lane (Cot) ecm enc) Gira}, 


The set Hf’, of all opportune delivery periods for a particular route, mission period, and 
cargo type combination 1s: 


ji”, = (P, | (P, x (7. 2h} © HW, {2 © P,}. 


In addition to the conditions mentioned above above, ib ee , and Vv may contain 
other implicit constraints which are specific to a given problem. The variables within the 
oo are restricted by sets of allowed combinations of primarv indices. The variables 

Seeecre delined over the derived set ’ , While the variables Y_., are defined over the 


derived set I’. The problem 1s formulated as a linear program: 


Minimize we Car Lines as Dye ae, ae ae » ea a 


(a,r,t,hye V ee Oeargh (1 H’ (r, z,h)e L’ 

Subject To 

eta RCAR HN SOASUP,, @HiES GL) 
(r,Aje a 

eS >, ent a Von + Nip = MREO,, (rn2A)eL’ (2.2) 

aé Ve Te W rth Te ie 

I4~7 44 
linear side constraints involving X",,,, and Y,,-7 (23) 
Ue ae (2.4) 


aie cS) 
0 <S iv <h; (2.6) 


The interpretation of the constraints in the Movement Requirement Reallocation 
Model is similar to that of the constraints in the Generalized Transportation Formu- 
lation. Constraint (2.1) ensures all movement requirements are satisfied by mission de- 
livery variables, using reallocation variables to transfer movement requirements into or 
out of the node which is represented by the constraint. Constraint (2.2) ensures the 
aircraft pool is not exhausted. Figure 2 shows a simple depiction of the network relation 
obtainable from these two constraint sets. The mission delivery arcs are shown with 
solid lines while the reallocation arcs are shown with dashed lines. Equation (2.3) re- 
presents the linear configuration constraints while equations (2.4), (2.5), and (2.6) are 
nonnegativity constraints. 

Since the network involved with the transformed model is no longer a bipartite 
graph directed from Sto L’, there is the distinct possibility of illegitimate flows to exist. 
An illegitimate flow exists when nonzero reallocations occur on a directed path of real- 
location arcs connecting three or more mission requirement nodes. In the event the 
ending node is the same as the originating node for the path, a cycle is created. A di- 
rected path implies that a node on the path will have both an arc importing flow into it 
and another arc exporting flow from it. If an illegitimate flow exists in the model, then 
a movement requirement is shipped with a lower total cost by using successive reallo- 
cations across intermediate time periods. This avoids incurring the larger cost of real- 
locating a movement requirement directly from the opportune delivery period to the 
nussion period. Since the resulting shipment by either reallocation process is the same, 
the use of the intermediate transfers violates the intended relative cost structure of the 
model. 

To prevent the occurrence of illegitimate flows in the model, a restriction on the 
costs for the reallocation variables and the infeasibilitv variables of the objective func- 
tion 1s required. Letting the subset { k,/,m } < P, represent any three time periods of the 
set of opportune delivery periods with both rv and fA being fixed, the following restriction 
will be imposed on the costs involving a movement requirement reallocation between 


periods & and o1: 


¥ y y 
Coo “Crap Cope (1) 
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Lemma: With regard to the Movement Requirement Reallocation Model, condition (1) 
is sufficrent to prevent the simultaneous tmportation and exportation of flow on reallo- 


Cation arcs at a L node 


Proof: Let node / be the L’ node with both route r and cargo tvpe A fixed and with the 
opportune delivery period of /. Let this node contain a directed path between nodes k 
and m. Condition (1) states the cost for using the infeasibility arc associated with node 
/ is less than the total cost of the two directed arcs connecting it with nodes A andm. 
Since the infeasibility arc associated with node / is unbounded and the Movement Re- 
allocation Model is a minimization problem, the infeastbility arc will always be used 
prior to the use of anv combination of two or more reallocation arcs through the di- 
rected path containing the node /. Therefore, the simultaneous importation and 


exportation of flow on reallocation arcs can not occur in the Movement Requirement 
Reallocation Model. QED. 


While condition (1) 1s sufficient to prevent illegitimate flows in the transformed 
model. additional conditions are required to ensure the use of variables in the models 


are Valid. The additronal conditions are: 
Cy a 2 
rkemh = ren (2) 
x y x a 
Cae < Cron te Ca ° (3) 


Condition (2) states the cost of using the reallocation variable for node k must be less 
than the cost of its infeasibility variable. This ensures reallocations of movement re- 
quirements occur before a movement requirement 1s determined to be undeliverable. It 
1s noted that when restrictions (1) and (2) are assembled together, a triangle inequality 
exists separated by the cost associated with the infeasibility variable. However, the tri- 
angle mequality 1s not sufficient by itself to prevent illegitimate flows since a directed 
path could exist between nodes &, /, m while no variable existed linking nodes & and m 
directly. The addition of the infeasibilitv cost to the restrictron prevents the stmultane- 
ous importation and exportation of flow for this type of directed path which would 
otherwise satisfy the triangle inequality. Restriction (3) 1s imphed by the definition of the 
opportune delivery period. For the definition to be meaningful, movement requirements 
must be satisfied within missions periods which are the same as the opportune delivery 


period prior to the reallocation of the requirement to a node with a different opportune 


delivery period. Consequently, the cost of satisfving the requirement in the same mission 
period as the opportune delivery period must be less than the cost of reallocating it to 
a different opportune delivery period and the associated cost of performing the mission 
during that period. 

The cost scheme in the General Transportation Formulation is convertible to an 
equivalent structure for the Movement Requirement Reallocation Model if cost re- 
striction (1) 1s satisfied for both models. It is assumed restrictions (2) and (3) already 
exist in a form useable to the General Transportation Formulation. Because the mission 
delivery costs in the models are independent of the reallocation costs, the usage penalty 
for aircraft are exclusively assessed on the mission delivery variables since they represent 
the actual shipment of cargo bv aircraft. Both the timeliness of delivery penalties and 
the cargo assignment penalties represent costs exclusively associated with the reallo- 


cation variables. Therefore, the transformed cost relations are: 


Cirh = APEN, 


Ch = EPPEN» + UEIPENe 


r 


When comparing the two models, it is important to note the distinct difference be- 
tween the sets L of the Generalized Transportation formulation and L’ of the Movement 
Requirement Reallocation Model. Both are used to restrict the first constraint set of 
their respective models. There will typically be many more constraints in the constraint 
set of equation (2.1) of the Movement Requirement Reallocation Model. This occurs 
because the index set L’ does not eliminate the node (een when aimovenient fe- 
quirement does not exist for those indices. This null movement requirement node must 
still exist since movement requirements from other opportune delivery periods mav be 
reallocated to it. This presents a tradeoff between the two models involving the re- 
duction of variables versus the reduction of constraints. It 1s expected there will be a 
cross-over point for the solvability of instances of problems dependent upon the density 
of the number of nonzero movement requirements in relation to the total possible 
number of movement requirements given the size of the primary index sets. The General 
Transportation Formulation should perform better for low density instances of prob- 
lems. The Movement Requirement Reallocation Model will improve its performance in 


comparison with the General Transportation Model. 
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D. INDEPENDENT TIME-WINDOW SEPARATION ALGORITHM 

For many operation plans involving the aircraft deployment problem, the size of an 
instance of the model may be intractable or have an unacceptably long solution time. 
This section of the thesis explores a specialized solution technique which can be used to 
solve a certain subset of these difficult problems. The subset of problems involves those 
cases where either the early deliver or the late delivery of cargo 1s unacceptable. When 
the model 1s restricted to allowing only earlv or only late deliveries of cargo, a simple 
algorithm can be used to find independent subproblems within the entire problem. 

Assuming a problem under consideration has no gross inadequacies in the available 
supply of aircraft with regard to the total amount of cargo to be shipped, it 1s anticipated 
that independent subgroupings of periods may exist. These subgroupings or time- 
windows form a subproblem which is independent of the remainder of the model if the 
movement requirements within the tume-window can be satisfied entirely using the air- 
craft available in that same time-window. It is usually difficult to tell where these in- 
dependent time-windows would occur for a particular instance of the problem without 
solving the entire instance. However, it is possible to iteratively define and solve sub- 
problems of independent time-windows when either earliness or lateness 1s not allowed. 

Consider an instance of the deplovment problem where early deliveries are disal- 
lowed. A time-window 1s independent of future periods if no late deliveries occur outside 
of this ttme-window. This happens in general because the cost structure of the deplov- 
ment models enforce the delivery of movement requirements as close to the opportune 
delivery period as possible. Since movement requirements can not be delivered in peri- 
ods earlier than the ttme-window and all deliveries to time periods later the time-window 
have higher costs, then all deliveries are made within the time-window if there are sufh- 
cient assets available. 

An iterative algorithm can be used to solve the entire problem which makes use of 
this property. The algorithm starts bv initially solving a subproblem containing a time- 
window with only the first period of the ttme horizon. A time-window of length 7 is 
w={t, f+1, 1+2,..., tt+n—1}, where r+heP and k=90. 1,...,n—1. If all 
deliveries are made within w , the subproblem is labeled as independent, the solution for 
the subproblem 1s saved, and q 1s revised to include only the next period of the time 
horizon. If not, the solution 1s disregarded and the next period of the time horizon is 
added to the existing period in w The next subproblem is solved based on the revised 
w. This process of solving the subproblem and then updating w continues until the last 


period of the tume horizon has been solved. Thus, a forward directed algorithm 1s ob- 


24 


tained which uses an ascending ordering of the periods (earliest to latest) to solve the 
entire problem by progressively solving subproblems over w . A similar procedure is 
possible when late deliveries are disallowed which iteratively solves subproblems based 
on a descending progression (latest to earliest) through the periods of the time-horizon. 
A simple algorithm outlining the general procedure is shown in Figure 3. 

This algorithm is particularly well suited to instances of the deplovment problem 
Where the demand for movement requirements fluctuates substantially from period to 
period or is unimodal in nature with respect to time. As the size of the problems at- 
tempted increases, the benefit in solving the smaller independent time-windows of those 
problems increases. Two other important benefits of the algorithm are its capability to 
model extended time horizons for certain instances, and its ability to be restarted from 
advanced positions in the time horizon. The worst case instance for this algorithm 
would have deliveries outside of the opportune delivery periods for all periods in the 
time-horizon. In this case, the algorithm will ultimately solve the entire problem, being 
no worse off for using the algorithm than having pointlessly expended effort solving the 
smaller subproblems. If the problem ts still intractable after using this algorithm, there 
is still valuable partial information which was obtained from iterating through the vari- 
ous time-windows of the algorithm. Even if no independent time-windows Were found, 
bottlenecks in the problem would most likely be identifiable from the partial information 
or an advanced starting position found for continuing the algorithm following trme- 
windows which are nearly independent. Another use for the algorithm is as a tool to 
determine bounds on instances of problems where both early and late deliveries are al- 
lowed. The procedure would be to solve the problem twice, once where earliness was 
disallowed and once where lateness was disallowed. Since either condition is a restriction 
of the original problem, either offers a bound on the unrestricted problem. By studving 


both solutions, valuable information for the total solution can be obtained. 


E. INTEGER ROUNDING MODEL 

In order to make the solutions from the previous models compauble with the exist- 
ing post-solution processor in the DAP, an integer rounding method 1s proposed. From 
the solution of any of the previous models, the amount of each movement requirement 
to be transported on the various aircraft tvpes over a specific route and mission period 
are determined. The number of aircraft required for each mission can be obtained from 


this solution using either of the following equations: 


tJ 
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ALGORITHM: Independent Time-Window Separator 


iVvEU Military Aircraft Deployment Problem formulated 
using a Generalized Transportation Formulation 
or Movement Requirement Reallocation Model. 
Periods P will be preordered from earliest to latest 
(ascending order). 

OUTPUT: List of aircraft mission assignments to 
movement requirement. 


{ Initialize: 
IF (no lateness allowed) 
reorder P in descending order 
w = time-window of the first period of the ordered P 
WHILE ( w e P) 
{ WHILE ( MREQw) = 0) w = next period 
Pec cs Fo 
{ SOLVE the subproblem over w 
Piw) = number of deliveries outside of w 
TEC Dias 0) 
{w = w + next period } 
ELSE 
| { w = next period } 


j 
OUTPUT 





Figure 3. Independent Time Window Separation Algorithm 
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where the variables on the right-hand side are at their optimal values. The value will 
often be fractional for AREQ,, . 

The rounding model of this study provides one method to round the continuous 
solution. This model attempts to induce integer solutions by rounding up fractional 
values of AREQ.,, provided there are unused aircraft available of the same aircraft type 
as the variable being rounded. This avoids leaving cargo undelivered when there are 
unassigned aircraft which could deliver it. The underlving principle is to sort the frac- 


tional part of the continuous solution by magnitude from highest to lowest for each 
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aircraft tvpe during a mission period. The fractional parts are then rounded up to the 
value of one, maintaining the ordering of the sort, until supphes for an aircraft tvpe 
during a mission period are exhausted. Any fractional part of a value which cannot be 
rounded up because of aircraft supply restrictions 1s onutted (amphcitly truncated to 
zero) from the solution. This principle is programmable using a variety of methods. 
The Integer Rounding Model which follows applies this principle in a linear programming 


environment. The primarv index sets of the Integer Rounding Model are: 


aircraft types. 
routes. 


period of aircraft missions. 


sie Bs 


Caneo types. 


The given data for the problem is: 


ASUP,, number of aircraft type a available to perform missions during period ¢. 


AREQ,,, continuous solution from previous model expressed as the number of air- 


craft used. 


i. fractional part of AREQ,,, 
The variables for the problem are: 
Z. indicator for rounding up a fractional aircraft of type a used on route r for 


a mussion during period /. 


In addition to the primary sets the following derived sets are necessary. The set S 
of all possible combinations of aircraft and mission periods such that aircraft are avail- 
able is the same as described in the previous models. The set f of the aircraft tvpes, 


routes, and mission periods such that fractional aircraft solutions exists 1s: 


f= {AOR xan, 2 0}. 


The problem is formulated as a linear program: 


Mi axinuze yy Lee Lae 


(a,r,Hef 
Subject To 
» Zart Ss ASUP,, = »y (AREQar, — Ears) (ae S (3.1) 
reR reR 
nS Hort) (3.2) 


The objective function maximizes the total of the fractional parts of the continuous 
aircraft solution Which are rounded up using the binary indicator variable Z,,,. Con- 
Straint (3.1) ensures the indicator variable Z,,, does not exceed the remaining supply of 
aircraft. The remaining supply of aircraft is shown in the right hand side of the equation 
as the total supplv of available aircraft less the integer portion of the continuous aircraft 
solution previously calculated. Constraint (3.2) 1s the zero-one restriction on Z,,, . 

The Integer Rounding Model, when solved as a continuous linear program replacing 
equation (3.2) with O< Z,,,< 1 , will always generate integer solutions. To understand 
this from a mathematical programming perspective, the nature of the constraints of 
equation (3.1) must be explored. Letting U,,, be the unused aircraft represented ss ame 


right hand side of equation (3.1), consider the matrix form of the relaxed integer model: 


max { EZ: AZeU Pee 


Where A 1s the technical coefficient matrix. Each of the constraints generated from 
equation (3.1) forms a set of ones for each row of the A matrix. Since the indices of a 
ands are fixed for each equation, mutually exclusive sets of ones are generated for each 
constraint across the routes. Therefore, there will only be a single one in anv column 
of A, the rest of its elements being zero. Each of the rows of the A matrix are part of a 
set of generalized upper bound (GUB) constraints. Additionally, a simple upper bound 
of one 1s applied to each variable. The proof of integralhity follows directly from the 
unimodularity of the A matrix and the integrality of both the variable bounds and U, 
which imphes any extreme point is integer [Ref. 13]. 

At this point, it must be emphasized that this model does not necessarily provide 


an optimal conversion of the original deployment model solutions to integer. It is an 


approximation. However, the approximation should be close to the optimal solution 
when nontrivial numbers of aircraft are required for most routes during a military oper- 
ation. Notice this condition is the same as that previously mentioned for the deployment 
models. Consequently, a reasonably exact measure using a rounded solution should be 


obtained when an instance for one of the deplovment model satisfies this condition. 
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Il. COMPUTATIONAL RESULTS 


This chapter describes the results of testing representative data on the models and 
the algorithm previously described. The data for the tests were obtained from 
ULSTRANSCOM. The particular data set is similar to what would be expected to be 
solved on the DAP. The models were tested using this data on a personal computer with 
similar characteristics as those which the DAP is currently implemented on. The com- 
puter used was an IBM Personal System, 2 Model 50 with 640k of addressable memory 
and a 10 Mhz 80287 math coprocessor. 

The data set involves four aircraft types, nine routes, twelve time periods, and two 
cargo types. The nine routes are spread over three embarkation ports on the east coast 
of the United States and three debarkation ports in Europe. The number of aircraft 
available during any period ranges from 186 to 354. The test problem requires the de- 
liverv of 83 different movement requirements involving 219,693 short tons of freight. A 
crude approximation for the number of aircraft trips required for this problem vields 
over 4000 round trips between the different ports of embarkation and debarkation. This 
large number makes the use of a continuous variable deployment model appropriate. 
For this particular data set, the maximum amount of lateness for any movement re- 
quirement is three days, while the maximum amount of earliness is only restricted to 
periods within the time-horizon. 

Considering the number of index combinations for routes, time periods, and cargo 
tvpes, there are 216 different possible movement requirements which could occur for this 
problem. The $3 movement requirements which actually exist out of the possible 216 
were considered to be representative of the relative density of the actual movement re- 
quirements to the total possible movement requirements for a problem. To examine the 
performance of the General Transportation Formulation (GTF) versus the Movement 
Requirement Reallocation model (MRR), the original set of movement requirements 
(data set ORIG) was expanded or reduced to yield a higher and lower density test in- 
stance. The higher density instance (data set HIGH) contained 135 movement require- 
ments involving 225,103 short tons of freight. The lower density instance (data set 
LOW) contained 33 movement requirements involving 106,906 short tons of freight. 
All other data elements remained as in the original data set. In order to test the Inde- 


pendent Time-Window Separation Algorithm (ITWSA), a second comparison run was 
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performed with lateness disallowed in the model. The ITWSA was applied to the MRR 
model. A test without the separation algorithm was also run for comparison. 

The solution characteristics for the linear programs of the various tests are shown 
in Table 1, and the total GAMS execution times and the linear program solver times are 
shown in Table 2 and Table 3. Test type number one is for test runs where lateness was 
allowed in the instance and test type number tio is for test runs where lateness was 


disallowed. The type of solution technique is indicated by the name of the model. 


Table 1. LINEAR PROGRAM SOLUTION CHARACTERISTICS 


MODEL DATA SET LOW | ORIG HIGH 


GTF-1 

equations 

vaniables 

nonzero coefficients 


MRR-1 
equations 
variables 
nonzero coefficients 


MRR-2 
equations 

vaniables 

nonzero coefficients 


ee S\-2 
equations 
variables 
meonzero coeflicients 
max eqns. iteration 
max vars. iteration 
miax nzcoef. iteration 


SeolATION: 



























GTF-1 Generalized Transportation Formulation -run | 
MRR-1 Movement Requirement Reallocation Model - run 1 
NIRR-2 Movement Requirement Reallocation Model - run 2 


ITWSA-2 Indep. Time-Window Separation Algorithm - run 2 





Table 2. TOTAL EXECUTION TIME (minutes) 


a [eee oe | 
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Table 3. LINEAR PROGRAM SOLVE TIME (minutes) 
MODEL DATASEY LOW | ORIG | HIGH | 
SGTFA 500 








2) 
~~] 
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The GTF required a consistently lower number of equations in comparison to that 
of the MRR model. This allowed the GTF to perform well for the lower density in- 
stance; however. as the densitv increased, its size and solution time increased sharply. 
The increase was so severe, that the GAMS program generator was unable to create the 
linear program for the higher density instance because of inadequate memory capacity 
in the computer. This sharp increase in problem size in relation to the number of 
movement requirements demonstrates that this model’s ability to solve the deployment 
problem ts highly dependent on the density of the movement requirements. This prop- 
erty 1s not desirable since a variety of instances of the deployment problem mav need to 
be solved. 

The MRR model had a shghtly larger size for the lower density instance than that 
of the GIF. The growth in model size was very moderate from the lower density in- 
stance to the original and higher densitv instances Which 1s in sharp contrast to the rapid 
growth found in the GIF. Despite the large disparity in the number of movement re- 
quirements used in these different instances, the MRR model remained stable in both 
size and solution time. This stability is directly related to the constant number of 
equations found in all instances of the problem, as well as the smaller quantity of vari- 
ables and nonzero coeflicients required when comparcdsto me onli. 

The TT WS.\ performed quite well for this particular data set. The algorithm was 
able to identifv eleven independent time windows in the twelve period problem. The 
accumulated linear program solution times were extremely low. The total execution time 
was stable: however, it was also larger than the nonseparated MRR model solution time. 
While the solution time was small, the total execution time was larger because of the 
additional overhead of generating all of the subproblems. This data set demonstrates the 


usefulness of the algorithm: however, it 1s not appropriate to think all data sets would 


perform nearly as well. The algorithm’s performance is based on the ability to find in- 
dependent trme-windows within a problem. For problems with few or no independent 
tume-windows, the solution time would be expected to increase substantially. Although, 
the maximum problem size would never be any larger than that which would be required 
to solve a problem without the separation algorithm. 

The Integer Rounding Model was used to transform the continuous solutions of 
previous model test runs to integer values. The rounding model performed well giving 
deployment plans which were consistent with the objectives indicated in the the data set. 
The quality results are directly related to the large number of aircraft required across the 
different routes and time periods. Additionally, it is noted that the time required to run 
the Integer Rounding Model was modest. 

The utility or supportability of the operation plan from the perspective of air trans- 
portation is indicated by the number of aircraft required and the amount of early, late, 
and undelivered cargo. For this particular data set, 14,247 stons of cargo were delivered 
one period early and 2,284 stons of cargo were delivered one period late. All of the re- 
maining 203.162 stons of cargo were delivered during their opportune delivery periods. 
If the periods are considered to be of one dav duration. then the shipment of all move- 
ment requirements used 1455 aircraft-days for the continuous variable solution and 1477 
aircraft-davs for the rounded solution. The overall utilitv of the operation plan would 
depend on these numbers being acceptable for the nulitary actions being considered. 

For comparison, approximations for this deployment problem were calculated using 
the ITWSA. When lateness was disallowed, 14247 stons of cargo were delivered one 
period early while 2284 stons went undelivered. When earliness was disallowed. 11195 
stons of cargo were delivered one period late and 5336 stons were delivered two periods 
late. There was no undelivered cargo when earliness was disallowed. Bounds on the 
amount of earlv, late, and undelivered cargo are easily deternuned from these solutions. 
Comparisons of the number of aircraft-davs required 1s not as straight-forward. As the 
amount of undelivered cargo increases the number of aircraft-days required decreases, 
but, as the amount of earliness or lateness increases, the number of aircraft-davs in- 
creases or decreases depending of the size of the aircraft and the priority of usage for the 
different aircraft types. For the current penalty scheme, the number of aircraft-days re- 


quired decreased for both the no-earliness and no-lateness approXimations. 
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IV. SUMMARY 


This thesis has presented techniques which can be applied to the aircraft deployment 
problem. A formulation based on the generalized transportation problem laid the 
foundation for exploring other methods. This formulation, in itself, was shown to be 
potentially unstable in size and solution time within the limitations of the solver and 
computer when solved as a linear program. This instability depended directly on the 
number of movement requirements in a given instance of the problem. Despite this, the 
model may still be useful with other optimization techniques such as decomposition, 
dynamic column generation. or network algorithms. The GAMS modelling program for 
this Generalized Transportation Formulation is given in Appendix A. 

A transformation of the Generalized Transportation Formulation was performed 
by separating the original transportation variables into those which were used for satis- 
{ving movement requirements in the time they were required and those which were used 
to satisfy requirements which were early or late. These variables were aggregated to 
produce an efficient reformulation of the former model. This model was found to be 
quite stable and performed adnurably on a personal computer. It is the model of choice 
when the problem 1s solved as a standard linear program. [The GAMS modelling pro- 
gram for the Movement Requirement Reallocation model is presented in Appendix B. 

A time window separation algorithm was proposed as an instrument to attack large 
problem instances. The model provides optimal solutions for problems which allow only 
early or only late deliveries of movement requirements. The model can also be used to 
find bounds on the solution for models allowing both earliness and lateness bv restricting 
the earliness or lateness variables. The model attempts to find independent groups of 
time periods within a problem so that the entire horizon of time periods does not have 
to be solved in its entirety for any iteration of the algorithm. Rather, a collection of 
subproblems are solved giving a set of solutions which can be assembled into the total 
problem solution. The MRR model in Appendix B is equipped with the Independent 
Time-Window Separation Algorithm. 

Lastly, an integer rounding model was proposed which allocates the fractional part 
of the continuous solution to the residual supply of available aircraft. This model 1s an 
effective tool to convert continuous variable solutions to integer when there are not 


substantial numbers of low values to be rounded. Although this model has integer 
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variables, 1t can be solved as a continuous linear program with guaranteed integrality. 
The model provided quality conversions of the continuous solutions of the former 
models for the data sets tested. The GAMS modelling program for the submodule of 


the Integer Rounding Model 1s presented in Appendix C. 


ee) 
ts" 


APPENDIX A. GAMS PROGRAM FOR THE GENERAL 
TRANSFORMATION FORMULATION 


Sore GENERAL TRANSFORMATION FORMULATION VERSION (6922. 


SONTEXT 
MICHAEL C. PUNTENNEY 


ACKNOWLEDGEMENT: THIS MODEL IS BASED ON THE ORIGINAL AIR 
LIFT ALLOCATION PROBLEM DESCRIBED IN 
THE TRANSCOM LIFT OPTIMIZER VERSION 
(88.6.14) BY RICHARD E. ROSENTHAL. 
SOFFTEXT 


SOPFUPPER  OFFSYMARER OFPSYMiism 


* DATA 
Neseeeen em eee e eee 
SETS 

Pm all periods in time horizon i. / 
A aircraft types y / 
H cargo types / | 
i ports ii i 
R routes / /; 


ALIAS( Pm, Po); 


SH eS: 
OCR,I) mapping of route to embarkation port 


D(R,I) mapping of route to debarkation port 
COM(A,H) compatible aircraft cargo pairings 
/ /; 


TABLE ASUP( Pm, A) quantity of available aircraft ; 

TABLE MREQ(R,Pm,H) movement requirement ; 

TABLE ALPCR em.) ordinate of available to load period for Pm ; 
TABLE RCAP(A,R,Pm) freight capacity of aircraft types on routes ; 


PARAMETER PCAP(I,Pm) freight throughput limit at ports ; 


SCALARS 
MAXLTE maximum allowable number of late delivery periods; 
MAXERL maximum allowable number of early delivery periods; 
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Omewaxanaeeeeeweweeeweesw eee 2 = = 


ee ee 


meevise KCAP and ASUP if zero aircraft of route capacity are available 
RCAP(A,R,Pm)SC(ASUP(Pm,A) EQ 0) = O-; 
ASUP( Pm, A)S$(SUM(R, RCAP(A,R,Pm)) EQ 0) = 


SET Q(A,R,Pm,Po,H) derived index set for allowed X variables; 
QCA,R,Pm,Po,H)SCOM(A,H) = YESS$(( ORD(Pm) - ORD(Po) LE MAXLTE ) 
AND ( ORD(Po) - ORD( Pm) LE MAXERL ) 
AND ( ORD( Pm) GE ALP(R,Po,H) ) 
AND ( MREQ(R,Po,H) GT 0 ) 
AND -(C ROAFCA, Ry em) Gl OF jee. 


ov mweme@eeFr@eeeese sea eweee2e 2] |] | = 


eo 
vm@awmw@e@ewese see eee ese es Ss es = 


PARAMETERS 
APEN(A) aircraft priorities 
HPENCH) EedrceO prtoOri ties 
TPEN( Pm, Po) time of delivery penalty from Po to Pm 
NPEN aircraft infeasibility override penalty ; 
SCALAR 
ASCALE scaling parameter for APEN 
HSCALE scaling parameter for HPEN 
SCALE Scaling factor £Or solution staba lity FAOORUO/ 


Default HPEN(H) assumes cargos entered in ascending priority: 
HPEN(H) = ORD(H) ; 


* Note: APEN(A) values specify order in which assets should be used, 
* di.e., highest value signifies most important to conserve. 
Default APEN(A) assumes assets entered in descending priority. 
APEN(A) = CARD(CA) - ORD(A) + 1 ; 


* Default TPEN(Pm,Po) 
TPEN(Pm,Po) $ ( ORD(Pm) GT ORD(Po) ) = 
2 + (ORD(Pm) - ORD(Po)) / CARD(Pm) ; 
TPEN(Pm,Po) $ ( ORD(Pm) LT ORD(Po) ) = 
2 + (ORD(Po) - ORD(Pm)) / sai - 1/( 2 * CARD(Pm)): 
TPEN(Pm,Po) $ ( ORD(Pm) EQ ORD(Po) ) = 


Determine penalty scale magnitudes ; 
HSCALE = ( 0.9 * 1/(2*CARD(Pm)) ) / CARD(CH) ; 
ASCALE = 0.9 * HSCALE / CARD(A) ; 


* Scale APEN and HPEN to fit penalty hierarchy 


ay! 


APEN(A) 
HPEN(H) 


ASCALE * APEN(CA) ; 
HSCALE * HPENCH) ; 


* Default NPEN greater than any penalty associated with Y and less 
than twice the minimum penalty associated with Y. 
NPEN = 33.57% 


Ea 
e@ee@e@eeqe@eeeeeeeeqedse2se2ee & 


a@e@eq@qoeedeeeoeeeeoere2e2 @ @ @= = 


VARIABLES 
X(A,R,Pm,Po,H) movement requirement assigned to aircraft mission 
NCR, Pm, i) infeasible movement requirements 
COST total penalty costs ; 


POSITIVE VARIABLES X,N; 


EQUATIONS 
BALANCE(R,Po,H) maintain freight balance constraints 
ASSETMAXCA, Pm) observe maximum available assets 
TROP Cie) observe throughput limits at ports 
ORIREE define objective function ; 
OBIDET a. 


SUM((A,R,Pm, Po,H)$Q(A,R,Pm,Po,H), X(A,R,Pm,Po,H) 
* (APEN(A) + TPEN(Pm,Po) + HPEN(H)S$(ORD(Pm) NE ORD(Po)))) 
+ SUM((R,Po,H), (NPEN + HPEN(H)) * N(R,Po,H)SMREQ(R,Po,H)) 
=F= COST - 


BALANCE(R,Po,H)SMREQ(R,Po,H) .. 
SUM((A,Pm)$Q(A,R,Pm,Po,H), X(A,R,Pm,Po,H)) 
+ N(R,Po,H) § MREQ(R,Po,H) 
=G= MREQ(R,Po,H) / SCALE ; 


ASSETMAX(A,Pm) S$ ASUP( Pm, A) 
SUNC CR Pe.) SOCASR 2 menos 
X(A,R, Pm, Po, Hy RGAPG emp eae CALE 
=L= ASUP(Pm,A) ; 


THRUPUT(I,Pm) $ ( PCAP(I,Pm) LT +INF ) .. 
SUMCRSO(R,I), SUM((A,H,Po)$Q(A,R,Pm,Po,H), X(A,R,Pm,Po,H)) 
+ SUM(RSD(R,I), SUM((A,H,Po)SQ(A,R,Pm,Po,H), X(A,R,Pm,Po,H)) 
=I= PCAP(i ba / SCAla.. 


) 
) 


Qo w@w@eweeeerwreedeoerveees wee @2@eoeea = 


MODEL SLIEIOR? y ALLy; 
OPTION LIMROW=0, LIMCOL=0, SOLPRINT=OFF, ITERLIM=5000; 
SOLVE LIFTOPT USING LP MINTMAING COst Se. 
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PARAMETER XACR,Pm,H,A) display aircraft mission cargo assignment 
PGR Seite OH) display late or early movement requirements 
PLANES(R,Pm,A) display aircraft assignment by mission; 

RACK yEM,H.,4) = SUMCPO, X. LCA,R,Pm,Po,H)) * SCALE ; 


XD(R,Pm,Po,H)$(ORD(Pm) NE ORD(Po)) = 
SUM(A, X. LCA,R,Pm,Po,H)) * SCALE ; 


N. LCR,Pm,H) = N. L(R,Pm,H) * SCALE ; 


PLANES(R,Pm,A) $ RCAP(A,R,Pm) = 
SUMCeHoCONCAL IH) Anh, bmeHeA) Je, RCAPCARR Pm) 5 
ORTION XA: 2: 3: 1; 
MEPEON XD: 2: 3: 1; 
@eetON PLANES: 2: 2: 1; 


DISPLAY XA, XD, N.L, PLANES; 


APPENDIX B. GAMS PROGRAM FOR THE MOVEMENT 
REQUIREMENT REALLOCATION MODEL 


USING THE INDEPENDENT TIME WINDOW SEPARATION ALGORITHM 


STITLE MOVEMENT REQUIREMENT REALLOCATION TRANFORMATION (IT) 
os WITH INDEPENDENT TIME WINDOW SEPARATION ALGORITHM 
i ( MAINFRAME VERSION - LOOPING SOLVE TECHNIQUE ) 
* VERSION (een 72c). MODUETED (eg. 2 7) 
SONTEXT 
MICHAEL C. PUNTENNEY 


ACKNOWLEDGEMENT: THE MRR TRANSFORMATION IS BASED 
ON THE ORIGINAL AIR LIFT ALLOCATION 
PROBLEM DESCRIBED IN THE TRANSCOM 
LIFT OPTIMIZER VERSIONSG@SSaG le 
BY RICHARD E. ROSENTHAL. THE ITWS 
ALGORITHM 1S BASED ON VERSION (87. 12) 2508 


SOFFTEXT 


SOFFUPPEKR OFESIMARER OPRS nro 


- 
ow @a&esese eeweewewertrervreeeeeweeeae = = 


DATA 
SEES 
Pm all periods in time horizon rs ve 
Tm( Pm) window of periods for current solve iteration 
TITER em) iteration counter ; 


ALIAS Eig EOnrE ear). 
ALIASC im, ie): 


SETS 
A aircraft types / y 
H cargo types / / 
i pOnts / / 
R routes / le 
SETS 


OCR,I) mapping of route to embarkation port 
/ i 
D(R,I) mapping of route to debarkation port 


COM(A,H) compatible aircraft cargo pairings 
/ . 


> 
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TABLE ASUP( Pm, A) quantity of available aircraft ; 

TABLE MREQ(R,Pm,H) movement requirement ; 

TABLE ALP(R,Pm,H) ordinate of available to load period for Pm ; 
TABLE RCAP(A,R,Pm) freight capacity of aircraft types on routes ; 


PARAMETER PCAP(I,Pm) freight throughput limit at ports ; 


SCALARS 
MAXLTE maximum allowable number of late delivery periods; 
MAXERL maximum allowable number of early delivery periods; 
DERIVED DATA 


Revise RCAP and ASUP if zero aircraft De route capacity are available 
RCAP(A,R,Pm)SCASUP(Pm,A) EQ 0) = O 
ASUP( Pm, A)S(SUM(R, RCAP(A,R,Pm)) EQ 0} = 0, 


SET WCR,Pm,Po,H) derived index set for allowed Y variables; 


WCR,Pm,Po,H) = YES $(( ORD(Pm) NE ORD(Po) 


( 
AND ( ORD(Pm) - ORD(Po) LE MAXLTE ) 
AND ( ORD(Po) - ORD(Pm) LE MAXERL ) 
AND ( ORD(Pm) GE ALP(R,Po,H)) 
AND ( MREQ(R,Po,H) GT 0 
AND ( SUMC ASCOM(A,H), RCAP(A,R,Pm) ) GT 0) ); 


PARAMETERS 
APEN(A) aircraft priorities 
Bee (CH) cargo priorities 
TEEN( Pm, Po) time of delivery penalty from Po to Pm 
NPEN aircraft infeasibility override penalty ; 
SCALAR 
ASCALE scaling parameter for APEN 
HSCALE scaling parameter for HPEN 
SCALE Sicalineetactor. fOr SOluLIon estab ierey J 1002007 
TOLER zero tolerance for defininig delay is Peo) yaa 


Default HPEN(H) assumes cargos entered in ascending priority: 
HPENCH) = ORDCH) ; 


Note: APEN(A) values specify order in which assets should be used, 
i.e. , highest value signifies most important to conserve. 

Default APEN(A) assumes assets entered in descending priority. 
APEN(A) = CARD(CA) - ORD(CA) + 1 ; 
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Default TPEN( Pm, Po) 
TPEN( Pm, Po) S33 GRD Pm) Cl 0RpGre a 
2 + (ORD( Pm) - ORD(Pe)) 7 CARDUEm), ; 
TPENCPm, Po)" S]@ ORDCPm) f1 ORE oe — 
2 + (ORD( Po) - ORD(Pm)) / CARD(Pm) - 1 / ( 2 *™ CARDC Pipes 


Determine penalty scale magnitudes ; 
HSCALE = ( 0.9 * 1/C2*CARD(Pm)) ) / CARD(H) ; 
ASCALE = 0.9 * HSCALE / CARD(A) ; 


Scale APEN and HPEN to fit penalty hierarchy 
APEN(A) = ASCALE * APEN(A) ; 
HPEN(CH) = HSCALE * HPEN(H) ; 


Default NPEN greater than any penalty associated with Y and less 
than twice the minimum penalty associated with Y. 
NPEN = 32553 


VARIABLES 
MAR > Pm cH) cargo assigned by aircraft mission 
YCR,Pm,Po,H) late or early delivry of movement requirements 
NCR,Pm,H) infeasible movement requirements 
COST total penalty costs ; 


POSITIVE VARIABLES A, Y, N ; 


EQUATIONS 
BALANCE(R, Pm,H) maintain freight balance constraints 
ASSETMAX(A, Pm) observe maximum available assets 
THRUPUT( I, Pm) observe throughput limits at ports 
OBJDEF define objective function ; 
ORIVEE a. 


aa (R,Tm,To,H), (HPENCH)+TPEN(Tm,To)) * YCR,Im,lo,H)SW(K, im) toe 
SUM( (R, Tm sH), (NPEN+tHRENCH) > SCR e in “H) SMREQ(R, Tne) 

: SUM((A,R Tm ,H)SCOM(A, H), APENCA) * X(A, lee, dive ,H)SRCAP(A, RK, ime 

=E= COST ; 


BALANCE(R,Pm,H)$Tm(Pm) .. 
SUM( ASCOM(A,H), XCA,R,Pm,H) $ RCAP(A,R,Pm) ) 
- SUM( Po $Tm(Po), Y(R,Pm,Po,H) $W(R,Pm,Po,H) ) 
+ SUM( Po $Tm(Po), Y(R,Po,Pm,H) S$W(R,Po,Pm,H) ) 
+ N(R,Pm,H)S$MREQ(R,Pm,H) 
=G= MREQ(R,Pm,H) / SCALE ; 


ASSETMAX(A,Tm)SASUP(Tm,A) . 
SUM( (R,H)SCOM(A,H), X(A,R,Tm,H)$RCAP(A,R,Tm) 
/ RCAP(A,R,Tm) ) * SCALE 
=L= ASUP(Tm,A) ; 


THRUPUT(1,Tm)$(PCAP(I,Tm) LT +INF ) .. 
SUM( RSO(R,I), SUM( (A,H)$COM(A,H), X(A,R,Tm,H)$RCAP(A,R,Tm) )) + 
SUM( RSD(R,I), SUM( (A,H)SCOM(A,H), X(A,R,Tm,H)SRCAP(A,R,Tm) )) 
=L= PCAP(I,Tm) / SCALE ; 


oeeweewewrtewrteeeewew @ ee @w w= = w= 
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MODEL LIFTOPT /ALL/ ; 
OPTION LIMROW=0, LIMCOL=0, SOLPRINT=OFF, ITERLIM=5000; 


PARAMETER XACR,Pm,H,A) display cargo assigned by aircraft mission 


OPTION XA: 


OPTION YA 


Sea 
ool 
wiitON NF: 2: 2: 1; 
SP IlON PLANES: 2 


WACK Sem, Po,H) display late of early movement requirements 
NE Chem yr) display infeasible movement requirements 
PLANES(R,Pm,A) display aircraft assignment by mission ; 


Ze ‘ 
2: . 


b 


‘ 


De aaa lp 


onmenweweerwwewuweeweeqeeeqeoeq®@ 2 & 


ogx wegeeqeeeqeeqeqeewFiFrFFw 2 = @&F Fs = = & 


PARAMETER 


ITWSA determines whether the Independent Time Window ; 
Separation Algorithm is ON or OFF 
( 1=0N; O=OFF ) 


ITWSA = 1; 


SET MREQTEST(C PM) terminate zero-MREQ iterations 
BLANKEND( PM) terminate last zero-MREQ solve 
SOLVIEST( PM) terminate solve loop when end test satisfied 


PALL( PM) all Tm if early and late deliveries allowed /T1*T12/; 
SCALARS 
PER iteration counter 
INFEAS determines infeasibilities for the current iteration 
CHECK test for nonzero movement requirement iteration 
TOTLCOST total cost of all feasible iterations / 0 / 
SIGN backward or forward iteration solve discipline 
SEND ending iteration count number ; 
SIGN = -1S(MAXLTE EQ 0) + 1S(MAXERL EQ 0) 
- 1$( MAXLTE EQ O AND MAXERL EQ O ); 
SIGN = SIGNSITWSA ; 
ITER = CARD(Pm)S(MAXLTE EQ 0) + 1S(MAXERL EQ 0) 
sts ( MAXLIE BOQ 0 AND HAXERG EO 0 ); 
ITER = ITERSITWSA ; 
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TEND = OS(MAXLTE EQ 0) + (CARD(Pm)+1)$(MAXERL EQ 0) 

- (CARD( Pm)+1)$(MAXLTE EQ 0 AND MAXERL EQ 0); 
Tm( Pm) SCSICNS ee) YES S(CORDC Page) CARDO Ems. 
Tm(Pm)S(SIGN EQ 0) PALL( Pm); 
Tm(Pm)$(SIGN EQ 1) YES $(ORD( Pm) EQ 1); 
TITER( Pm)$(SIGN EQ -1) YES $(ORD( Pm) EQ CARD(Pm)); 
TITER(Pm)$(SIGN EQ 1) YES $C(ORD(Pm) EQ 1); 


ou i tt al 


SOLVGES? (Pm) =) yi: 
LOOP TEESSOLVIEsi tLe). 


MREQTESTCPH) = YEos CCARDCD ete. 
BLANKEND( TITER) = YES ; 
LOOP (PESHRE@T ES TCE E) 


CHECK =  1S(SUM( GRO tm) 4 MR EOCR. wine ae 
WG Rais = ITER +rsiGns( Cubek NE oe, 
Tm( Pm) = YESSCORD Crm em Tek) a 


MREOTESMOCPH) S CHECK saga 
OR TPRRSEOS TEND) = NO, 


BLANKEND( TITER) S(SIGN*(ITER)— TEND) GE 0) s—siGe. 
LOOCP(TITERSELANREND( (1 ieee 
SOLVE LIFTOPT USING LP MINIMIZING COST ) ; 


MREOTEST CRMs — iis: 


ITER =" [ter + Slane. 
INFEAS = SUM((R,To,H), N.L(R,To,H)) ; 
TOTLCOST = TOTLCOST + COST. L$(TOLER GT INFEAS 

OR ITER EQ TEND) ; 
TITER(Pm) = YESS(ORD(Pm) EQ ITER ) ; 
Tm( Pm) = Tm(Pm)$(INFEAS GT TOLER) + TITER(Pm) ; 


SOLVIEST( PM)S(SIGN=#( 17ERS= TEND GE Cpe eee 


XACR,Pm,H,A) = X.L(A,R,Pm,H) ™* SCALE ; 
YA(R,Pm, Po,H) = Y.L(R,Pm,Po,H) * SCALE ; 
NF(R, Pm,H) = N.L(R,Pm,H) * SCALE ; 


PLANES(R, bm, 4) °S KCAP( ASRS) s) — 
SUM( HSCOM(A,H), XACR,Pm,H,A) )3/ SREAE GA Seine 


DISPLAY XA, YA, NF, PLANES, TOTLCOST; 


a 


APPENDIX C. GAMS PROGRAM FOR THE INTEGER ROUNDING 
MODEL 


SONTEXT 


INTEGER AIRCRAFT ROUNDING MODEL (88. 06.01) 


MICHAEL C. PUNTENNEY = MODTE TED CS9. 2, is) 


ACKNOWLEDGEMENT: THIS MODEL IS BASED ON THE ORIGINAL AIR 
LIFT ALLOCATION PROBLEM DESCRIBED IN 
THE TRANSCOM LIFT OPTIMIZER VERSION 
(3712523) BY RICHARD’ Eo ROSENTHAL. 


SOFFTEXT 


SPePPUPPER OFFSYMXREF OFFSYMLIST 


Wweeenreereececeee ee ew = = 


PARAMETER FRAC(CA,R,Pm) fractional PLANE values 
PLANCR, Pm, A) rounded plane values for deployment plan ; 


FRAC(A,R,Pm) = PLANES(R,Pm,A) - TRUNC(PLANES(R,Pm,A)) ; 


VARIABLES 
TOTAL total of fractional values rounded up 
PCA, R,Pm) rounded fractional values of PLANES ; 


Z.UP(A,R,Pm) = 1; 
Z.LO(A,R,Pm) = 0; 
7. FX(A,R,Pm)$(PLANES(R,Pm,A) EQ 0) = 


EQUATIONS 
ONE objective function 
ALLOCATE( Pm, A) allocate fractional values observing supplies ; 


OoJr  .. TOTAL =E= SUM((A,R,Pm), FRAC(A,R,Pm) * Z(A,R,Pm) 
$(FRAC(A, eel) NE Oy ; 


ALLOCATE(Pm,A)SASUP(Pm,A) .. SUM(R, Z(A,R,Pm)) =L= 
ASUP(Pm,A) - SUM(R, PLANES(R,Pm,A) - FRAC(A,R,Pm)) ; 
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MODEL ROUNDER /OBJF, ALLOCATE/ ; 
OPTION LIMROW=0, LIMCOL=0, SOLPRINT=OFF; 
SOLVE ROUNDER USING LP MAXIMIZING TOTAL; 


PLAN(CR,Pm,A) = TRUNCCPLANES(R,Pm,A)) + Z2.LCA,R,Pm) ; 


OPTIGN PLAN: 1. 2:35. 
DISPLAY PLAN ; 
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